Universal linear relations between susceptibility and T c in cuprates 
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We developed an experimental method for measuring the intrinsic susceptibility x of powder of 
cuprate superconductors in the zero field limit using a DC-magnetometer. The method is tested 
with lead spheres. Using this method we determine x f° r a number of cuprate families as a function 
of doping. A universal linear (and not proportionality) relation between T c and x ls found. We 
suggest possible explanations for this phenomenon. 



I. INTRODUCTION 

Among the basic properties of a superconductor is its 
ability to expel a magnetic field, i.e. the Meissner effect. 
In all the metallic superconductors the diamagnetic effect 
is complete, and below T c , the susceptibility, x, equals 
— f. In the cuprates high temperature superconductors 
(HTSC) the situation is far from being so simple, and 
there is growing evidence of samples showing incomplete 
Meissner effect and even paramagnetic Meissner effect 
. At the same time there is an accumulation of results 
showing that the superconducting ground state in these 
materials is inhomogeneous [2|. Therefore, it is possible 
that the partial Meissner effect (\ < —I) in the cuprates 
is an intrinsic property. This possibility motivated us to 
perform a comprehensive study of DC-susceptibility in 
cuprates. We look for correlations between T c and Xi m 
different HTSC families, and various doping. 

It is important to mention that Panagopolous et al. 
measured the AC susceptibility of La2-ySr y CuC>4 and 
HgBa2Cu04+,5 families 3]. However, they were not inter- 
ested in comparing the absolute value of x between the 
families and concentrated only on comparing the tem- 
perature dependence of the penetration depth between 
X and /iSR measurements, which resulted in very good 
agreement. 

The difficulty of determining the absolute value of x is 
caused by the granular nature of the cuprates, and their 
ability to pin flux very easily. Consequently, the mag- 
netization in these samples depends very much on the 
measurement procedure. For example, cooling a sam- 
ple in a field, or cooling in zero field and then applying 
the field, will result in a different magnetization. On the 
other hand, the intrinsic susceptibility of a sample must 
be well defined and one should be able to compare dif- 
ferent samples. 

Therefore, we first develop the condition under which 
the measurements lead to the intrinsic susceptibility of 
the cuprates. The development of these conditions is 
based on experience gained while trying to measure the 
magnetization of a bundle of lead spheres. Second, we 
look for correlation between T c and x, m different HTSC 
families with various doping. Our major finding is a uni- 
versal linear relation between T c and x- 

Our x measurements are done on a set of 
HTSC families, which are different in many senses. 



The different families are La2- 2/ Sr J/ Cu04 (LSCO), 
YBa2Cu30 y (YBCO) and its less known "cousin" 
(Ca !C Lai_ a; )(Bai.7 5 _ x Lao.25+x)Cu30j / (CLBLCO) sys- 
tem with 4 different values of x. The CLBLCO system 
in particular is ideal for our study due to several inter- 
esting properties. Each value of x generates the full su- 
perconductivity dome from the under-doped to the over- 
doped, and the maximum T c is x dependent. Thus, each 
x can be considered as a superconducting family. For all 
values of x and y CLBLCO is tetragonal, so there are 
no structural transformations that can cause a change in 
the volume of the unit cell. 



II. SAMPLE PREPARATION 

Ceramic samples of LSCO, YBCO and CLBLCO were 
made by solid-state reaction. For LSCO, stoichiomet- 
ric amounts of La203, SrC03 and CuO were mixed and 
ground using a ball mill. The mixtures were fired in air 
for 1-2 days; this was repeated three times. After pel- 
leting, the samples were sintered in O2 for about 64h. 
The sintering temperature varied between II00°C and 
1175°C, depending on the Sr level; then the samples were 
cooled to room temperature at a rate of 10°/h. 

For YBCO, the starting materials were Y2O3, BaC03 
and CuO. The mixture was fired in air at 910°C, then 
pelletized and fired again at 930°C; the last step was 
then repeated. We also prepared pellets using YBCO 
that was supplied by PRAXAIR. This sample is made 
by combustion spray pyrolysis; the grains' average size is 
3.9 //m. The pellets of the two kinds of YBCO were then 
sintered in 2 for 60h at 970°C, cooled at a rate of I0°/h 
down to 510 °C, and at a rate of 5°/h to 4I0°C. The 
samples were kept at 4I0°C for 5 days and then cooled 
down to room temperature at a rate of I0°/h. 

The results presented in this paper for YBCO are from 
the two types of samples. No difference can be detected, 
meaning that the results are not sensitive to the prepa- 
ration method of the samples. 

The samples oxygen level, y, was then reduced by bak- 
ing the sample in O2 and quenching the samples in liquid 
nitrogen. For very underdoped samples the reduction 
was done in nitrogen atmosphere. The reduction tem- 
peratures are listed in table [H] . The preparation of the 
CLBLCO samples is described elsewhere 



2 





y 




atmosphere 


Material 


92(l)K 


6.983 






PRAXAIR 


86.7(2)^ 


6.855 


530°C 


o 2 


Technion 


56.7(2)^ 


6.549 


740° C 


2 


Technion 


50.5(2)^ 


6.489 


8WC 


o 2 


Technion 


40(2)^ 


6.399 


840° C 


2 


PRAXAIR 


20(2)^ 


6.3 


580° C 


N 2 


PRAXAIR 



TABLE I: Summary of all the YBCO samples and the pa- 
rameter values used in their preparation. T r is the reduction 
temperature. 
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FIG. 1: (Color online) The phase Diagram of CLBLCO, 
LSCO, and YBCO after conversion of chemical doping to hole 
doping p using Eq. Q 



The oxygen level of all the samples was determined by 
iodomctric titration. In the LSCO samples the deviation 
of the Oxygen level from 4 is less than 0.005. 

The T c of all the samples is determined using resistiv- 
ity measurements. In order to compare all the samples, 
and to overcome the nontrivial problem of the relation 
between the chemical doping and p, we plot in Fig. ^ a 
unified phase diagram using the Presland et al. formula 
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T c /T c , max = 1-82.6( P -0.16) 2 



(1) 



which relates T c and the holes density p. 

Scanning electron microscopy (SEM) pictures show 
that the grain sizes in our YBCO, and CLBLCO for dif- 
ferent values of x, are of the same order of magnitude, 
and that the grains are agglomerates of crystalline whose 
typical length is I ~ 1 — 10 \im. An example can be seen 
in Fig. [3 These properties ensure that the demagneti- 
zation factor is similar for the different families. 
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2: SEM picture of a CLBLCO sample with x = 0.4 and 
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III. EXPERIMENTAL METHOD 

The susceptibility measurements were done using a 
home built magnetometer based on a primary coil, two 
compensating secondary coils, and an extraction motor. 
Some results were verified with QD-SQUID at Bar-Ilan 
University and with a Cryogenic S600 SQUID magne- 
tometer recently installed in our lab. The measurements 
were done in field cool conditions (FC), namely, for field 
changes the sample was warmed above T c and cooled 
down in the new field. Since we use a superconducting 
magnet there is always trapped flux in the magnet lead- 
ing to a constant shift in the field values. For that reason 
the magnetization is measured over a range of positive 
and negative fields. The susceptibility is defined by 



1 dm 

y = lim . 

A h^q V dH ' 



(2) 



where m is the magnetization obtained from the induced 
signal at the secondary coils, H is the external field, and 
V is obtained from mass/density. The calibration of the 
magnetometer is explained below. The definition of V 
requires clarification. The problem in powder samples is 
to achieve conditions where the volume out of which the 
field is expelled, V sc , equals V. The zero field cooling con- 
dition (ZFC) could result in a shielding volume V sc which 
is bigger than V because in certain geometries Josephson 
connections can lead to shielding currents enclosing non- 
superconducting regions in the sample. The field cool 
conditions, on the other hand, lead to the Meissner vol- 
ume and could result in a V sc which is smaller than V 
due to flux pinning. Therefore, our first challenge is to 
find the appropriate measurement conditions, where V 
obtained from mass/density is exactly the volume out of 
which the field is expelled for powders. 
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FIG. 3: (Color online) The volume fraction / of lead, which is 
the superconducting volume obtained from \ measurements, 
divided by the real volume taken from mass/density, plot- 
ted vs. the mass of samples. In the top axis we show also 
the dimensionality of the powder as the height it occupies 
in the sample container divided by its diameter. The solid 
squares represent sphere shaped grains and open circles rep- 
resent pancake-shaped grains. In the inset we show the mag- 
netization curves for three characteristic cases described in 
the text. 



In order to gain experience we performed a prelimi- 
nary experiment with Pb spheres where the theoretical 
X is well known since Pb is a type I superconductor, so 
A is negligible. We used sphere diameter of 0.5 mm, and 
assumed that \ is the susceptibility of a single sphere in- 
cluding the demagnetization factor (—3/2) and obtained 
V sc = lirriH_»o ~ j§ . We calibrated the susceptometer us- 
ing a few spheres mixed with sand so that they were very 
well separated from each other. Raw data are presented 
in the insert of Fig.[3]where we show curves of M — m/V 
vs H for 3 samples: (I) the few Pb spheres mixed with 
sand (17.24mg); (II) a layer of Pb spheres (55.81mg); and 
(III) a full container of Pb spheres (634. 5mg). In all cases 
a linear field dependence is observed at low fields. In the 
first and third cases we find the same slope at H — > 0, 
but both are different from the second case. This means 
that isolated spheres and a full container of spheres give 
the same result . 

Our findings in terms of / — V sc /V are summarized in 
Fig. where / is depicted as a function of sample mass, 
and as a function of height of spheres in the container (h) 
over its diameter (d), on the lower and upper abscissa, 
respectively. For a small number the Pb spheres, which 
form a 2D layer at the bottom of the container (h/d < 
1), we find / > 1. As the number of spheres increases, 
the volume they occupy in the container becomes 3D in 
nature (h/d > 1), and / converges to 1. We repeated 
the experiment with "pancake" shaped pieces of lead; 
the results are qualitatively the same. This leads to one 



FIG. 4: (Color online) Field cooled magnetization vs tem- 
perature in a variety of fields, in the small field limit where 
H <C $o/d 2 and d is typical grain size. 

of the findings of this work. As long as FC conditions 
prevail and we use large values of h/d, we can safely 
assume that V sc = V. This means that the magnetic 
field wanders inside the sample, in between the different 
grains, and fills all the empty spaces. 

In the cuprates x °f course is unknown, yet it is pos- 
sible to check if the experimental conditions developed 
for Pb apply here as well. For this we determined the 
magnetization m in an FC procedure for various ap- 
plied fields. We used fields which are small enough that 
even one flux quanta $0 = 20 Oe-/zm 2 cannot pene- 
trate our grains (cross section scale A ~ 1/j.m 2 ), namely, 
H <C $0/^ = 20 Oe. One such measurement is shown 
in Fig. Mostly data in the sub-Oe fields are com- 
pletely presented. The magnetization at the lowest tem- 
perature as a function of field is then plotted in Fig. [5] 
In this figure a single line seems to fit the entire field 
range. However, when zooming in on the sub-Oe region, 
which is shown in the inset, a global shift of the line with 
respect to the fit is seen between negative and positive 
fields (due to bias currents in the power supply). There- 
fore, we fit m vs. H to two different lines in a 10 Oe 
field range around zero magnetization, and obtain the 
susceptibility only from the averaged slope according to 
Eq. |21 However, outside this 20 Oe field range a kink in 
the magnetization appears which we believe indicates the 
first vortex that enters into a grain. Therefore, all data in 
our experiment were acquired using this 20 Oe field range 
in steps of 1 Oe. We are aware of works showing a signif- 
icant non-linear field dependence of the magnetization in 
single crystals, especially in very low fields (mOe) ja, 0] • 
However, we did not see any deviation from linearity in 
our experimental conditions over this 20 Oe range. 

To demonstrate that it is the intrinsic susceptibility 
of the cuprates that we are measuring, we present four 
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FIG. 5: (Color online) The zero temperature magnetization 
after a field cooled is plotted as a function of applied field. The 
field scale is shifted due to flux trapped in the superconducting 
magnet. Nevertheless, a straight line seems to fit the data 
well. Only a zoom in the zero magnetization region, depicted 
in the inset, shows that the shift is not identical on the two 
sides of zero magnetization. The susceptibility is determined 
by fitting the data to two lines, on different sides of zero 
magnetization, and taking the averaged slope. 



tests. First, we performed the susceptibility measure- 
ments as a function of mass for a CLBLCO sample with 
T c = 42.3 K in a cylindrical sample holder of 5 mm inner 
diameter. Again, large mass means a 3D cylindric like 
sample. In contrast the sample resembles a disk when 
the mass is small. As can be seen in Fig. HQ \ decreases 
with increasing mass and saturates. All our measure- 
ments are therefore done with large mass. Second, we 
use a set of sieves, and divide the powder grains into two 
groups: 20^m < d < 40/zm and d < 20/im where d is the 
characteristic size of a grain. We measured \ of these 
two samples both in FC and ZFC conditions, and the re- 
sults are shown in Fig. Ufa) and (b), respectively. There 
is hardly any grain size dependence in the FC measure- 
ments, especially when compared to the ZFC experiment. 
This indicates that the grain size does not play a role in 
determining x &s long as we use FC procedure. 

Third, due to the combination of weak flux pinning 
(compared to low T c superconductors) and high temper- 
atures, the time dependence of the magnetization can 
be very complex. Our main interest here is to find the 
optimal cooling scheme in a field in order to obtain repro- 
ducible magnetization at base temperatures. We checked 
the susceptibility of a sample as a function of the cooling 
rate. We found, in agreement with previous works [J 
that in FC conditions it is important to pass through T c 
slowly. Therefore, in all our measurements we cool the 
samples slowly enough so that no difference in the mea- 
surements is observed by cooling them even more slowly. 



FIG. 6: (Color online) \ vs. mass for a CLBLCO sample with 
T c = 42.3 K. 
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FIG. 7: (Color online) The magnetization curves for two grain 
sizes both in (a) FC and (b) ZFC. (c) depicts measurements 
for two different cooling rates. 

This is demonstrated in Fig. Etc) where we show that 
cooling at two different rates does not vary our result. 

As a final test we measured the magnetization of a 
T c ^40K LSCO sample, first in the form of a sintered 
pellet and then of the powder after pulverization. The 
results are shown in Fig. |H1 In the ZFC measurements 
there is a great difference between the magnetization of 
the pellet and of the powder. While for the powder we 
observed the linear behavior we saw before, for the pellet 
we find a more complex curve. Up to 20 G the calcu- 
lated susceptibility is almost -1, indicating a shielding 
supercurrent that keeps the entire volume of the sample 
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FIG. 8: (Color online) Magnetization measurements in FC 

and ZFC conditions of a sintered pellet and powder obtained FIG. 9: (Color online) The negative susceptibility as a func- 
by pulverizing the pellet. In FC conditions there is no differ- tion of hole doping p taken from Eq[U 
ence between the two samples. 



free of magnetic flux. Above this field the susceptibility 
decreases and reaches a value similar to that of the pow- 
der. The difference between the two samples is only the 
connection between the grains, those can be described 
as Josephson junctions with some average critical field, 
Hj c %. Above this field the inter-grain links can not sup- 
port the shielding current and flux penetrates into the 
space between the grains and we get local shielding of 
the grains as in the powder. 

On the other hand, the FC measurements give a differ- 
ent picture. The magnetization is linear in all fields, both 
for the pellet and for the powder. Furthermore, the sus- 
ceptibility is identical for both samples. This indicates 
that the inter-grain links cannot support any Meissner 
currents at all. The field is not expelled from the volume 
in between grains even at fields below Hj c %. 

The different behavior of the FC and ZFC measure- 
ments in the pellet sample demonstrate another advan- 
tage of our measurement procedure; it is not sensitive to 
the connectivity between grains. 

We interpret the results of all the above tests as reach- 
ing experimental conditions where small variations of 
these conditions have no affect on x- Therefore, we be- 
lieve that our experiments are in the limit where x is the 
intrinsic susceptibility of the cuprates. 



IV. RESULTS 

In Fig. we show the FC susceptibilty of all our sam- 
ples as a function of p, the hole concentration, where p is 
calculated using Eq. ^ The curves of — x vs. p resemble 
the phase diagram of Fig. 2] leaving no doubt that T c 
and x are somehow related. 



In Fig. HOf a) we present T c versus — x for all samples. 
We find that T c increases linearly (at low doping) with 
increasing — x, and the linear relation is identical for all 
families (within experimental errors). This is the main 
and theory-independent finding of this work. It is im- 
portant to mention that no correlation between x an d 
T c was found when the measurements were done in ZFC 
conditions. 



V. DISCUSSION 

The fact that all the samples obey the same linear 
relation between T c and x is very surprising, given the 
differences between these cuprate families. It may be 
that, because of the complexity of these materials, a new 
effective media theory is needed to explain this relation. 
Nevertheless, we would like to offer a simpler explanation 
for our data based on two experimental observations. On 
the one hand, it is well known that in the cuprates there 
is a universal linear relation between T c and the inverse 
in-plane penetration depth squared , known as the Ue- 
mura relation. This relation was revealed by a compre- 
hensive comparison of the penetration depth, measured 
by the muon spin relaxation (/iSR) technique, between 
different families of HTSC. On the other hand, based on 
the growing evidence for inhomogeneity in the cuprates 
and our observation that x is independent of grain size 
and connectivity (see Fig.|S|and|7J|, it is conceivable that 
the length scale of grain sizes observed in Figgis not the 
correct grain size. Therefore, we speculate that the ag- 
glomerates seen in Fig.|5]are made of a very large number 
of even smaller units stuck together, and that their num- 
ber is so large that the size of each one is smaller than the 
penetration length, at least in the low doping regime. In 
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FIG. 10: (Color online) (a) T c vs —% at T = 1.6 K for various 
samples of CLBLCO, YBCO and LSCO. (b) T c vs the muon 
depolarization rate a at T = 1.8 K for the same CLBLCO 
samples. Data for YBCO and LSCO are from Ref. Q 

this approach the true effective grain size length scale a 
would be a parameter to be determined experimentally. 

In type II superconductors such as the HTSC, the 
penetration length plays an important role in the sus- 
ceptibility, and x = — b(l — g(x)/x) 0, [n|, where 
x = a/X. In spherical, plane, and cylindrical shaped 
grains b = 3/2, 1, 1, and g(x) is the Langevin, hyperbolic 
tan gent , and the modified Bessel functions, respectively 
[Tfl lll| . Under the assumed long penetration length as- 
sumption (a/A < 1) g(x)jx can be expanded to give 
X = -ca 2 X~ 2 with c = 1/10, 1/3, 1/4 for the spheri- 
cal, plane, and cylindrical shaped cases, respectively. We 
further speculate that x(A) = 0(A~ 2 ) for all geometries. 
The anisotropy of the cuprates result in the replacement 
A = 1.3A a {, [T3 ■ After averaging over all grain shapes 
and sizes we expect 

_ a 2 

X = -C-ry (3) 

where c is the averaged c and is a number on the order 
of unity (including the factor 1.3). 

We also performed ^SR measurements on sintered pel- 
lets made from the same CLBLCO powders, at the Paul 
Scherrer Institute (PSI) Switzerland, by field cooling in 
3 kOe to 1.8 K. A full account of these measurements in 
CLBLCO is given in Ref. |I3, and the YBCO and LSCO 
data were taken from Ref. [14| . Figure HOf b) depicts T c 
versus a for all the samples. Here we used T c from /iSR 
as in the original Uemura plot. 

A comparison between the two plots reveals interesting 
information. First, by comparing the /iSR and suscepti- 
bility results we can estimate a. For this we fit both T c 
versus a and x m the underdoped region to straight lines 



with offsets <r g and x . The solid lines in Fig. HOf a) and 
(b) are given by 

T c = -K x ( X + Xo), (4) 

and 

T e = K^a + a ), (5) 

respectively, where = 62(5) K-/xsec and where K x = 
145(5) K. We determine a by making K^a and — K x x 
agree with each other once they are expressed in terms 
of A afc . Taking a = 7 x 10 6 A" h 2 [H where a is m \i sec 
and \ a b in Angstrom, and x from Eq. 0]we obtain 

7 x UfK^Xj = K x ca 2 X ab - 2 . (6) 

Solving this equation with c < 1 we find a > 200 nm. 
This length scale, which is smaller than the typical crys- 
talline size estimated from SEM, could be due to defects 
or an intrinsic separation into domains. The same length 
scale was also found independently by ac-susceptibility 
in YBCO and was ascribed to twinning [j^. However, 
our experiment shows that this is not the origin of a since 
CLBLCO and LSCO have no twinning. 

Second, there is an offset in both x an d o so that 
at X a b — * 00 we find T c ~ 10 K. This universal devi- 
ation from strict proportionality between T c and A~ 2 
is in agreement with the measurements of Zuev et al. 
on YBCO films 0|- The susceptibility offset could 
be explained by free spin that are present in under- 
doped HTSC and freeze as a spin glass 18] . The ex- 
pected susceptibility of paramagnetic spins is given by 
X = AirN /j, 2 j: j / (SksTV) , where the 4-7T is introduced here 
since we normalized the susceptibility in Fig. ^| so that 
X = — 1 for a superconductor [instead of — l/(47r)]. Tak- 
ing n e ff = 1.9/ib per Cu, T = 1.6 K, N = 3 spins in a 
unit cell, and V the volume of a cell we find x° ~ 0.1. 
However, free spin can not explain the offset in the fiSK 
a since they tend to increase a rather than decrease 
it, namely, with spins a is never zero. A different ex- 
planation for the offset, suggested in Ref. 01 j i s that 
T c oc A7i P with p ~ 1/2. This power law is most pro- 
nounced in the region T c < 10 K. This region is out of 
the scope of our measurements, but p ~ 1/2 at ultra low 
doping will give an artificial offset of the T c versus — x or 
a for "normal doping" (T c > 10). 

Third, panel (b) shows the well known boomerang ef- 
fect in YBCO and LSCO, namely, overdoped samples 
have higher a than underdoped ones with the same T c . 
In CLBLCO there is an anti-boomerang in both /iSR and 
X measurements especially for the x = 0.1 sample. This 
is a surprising result since it means that in overdoped 
CLBLCO, where the hole concentration is large, there is 
in fact a smaller superfluid density (n s cx A~ b 2 ) than in 
underdoped samples with smaller hole concentration. 
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VI. CONCLUSIONS 

We found a universal linear dependence for under- 
doped HTSC between T c and \ f° r different families, with 
doping as an implicit parameter. A possible explanation 
for this dependence is that in underdoped compounds the 
penetration depth A a {, is longer than an effective grain 
size length scale a, which is much smaller than the grains 
size measured using SEM. In that case \ is proportional 
to a 2 X~ b . By comparing T c (x) and T c (a) we estimate 
a > 200 nm. The amazing aspect of this new grain size 
is that it is independent of sample preparation, type of 
compound, and doping. It appears to be similar to do- 
main size in ferromagnets which are not determined by 
the sample size. In addition, our universal line does not 



cross the origin in the T c , \ plane indicating universal de- 
viation from strict proportionality between T c and A~ 2 . 
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